PROCEEDINGS 

OF SCIENCE 



o 

(N 

> 
O 

X 

6 



(N 
O 
(N 



Effects of Residue Background Events 
in Direct Detection Experiments on 
Determining Properties of Halo Dark Matter 



Chung-Lin Shan* 

Department of Physics, National Cheng Kung University, Tainan City 70101, Taiwan, R.O.C. 
Physics Division, National Center for Theoretical Sciences, Hsinchu City 30013, Taiwan, R.O.C. 



E-mail: clshan@mail . ncku . edu . tw 



We reexamine the model-independent data analysis methods for extracting properties of Weakly 
Interacting Massive Particles (WIMPs) by using data (measured recoil energies) from direct Dark 
Matter detection experiments directly and, as a more realistic study, consider a small fraction of 
residue background events, which pass all discrimination criteria and then mix with other real 
WIMP-induced signals in the analyzed data sets. In this talk, the effects of residue backgrounds 
on the determination of the mass of halo Dark Matter particle as well as on the reconstruction of 
its one-dimensional velocity distribution function will be discussed. 
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1. Introduction 

In our earlier work on the development of 
model-independent data analysis methods for 
extracting properties of Weakly Interacting Mas- 
sive Particles (WIMPs) by using measured re- 
coil energies from direct Dark Matter detection 
experiments directly [[|-|3]], it was assumed that 
the analyzed data sets are background-free, i.e., 
all events are WIMP signals. Active background 
discrimination techniques should make this con- 
dition possible. For example, by using the ratio 
of the ionization to recoil energy, the so-called 
"ionization yield", combined with the "phonon 
pulse timing parameter", the CDMS-II collab- 
oration claimed that electron recoil events can 
be rejected event-by-event with a misidentifi- 
cation fraction of < 1CT 6 . The CRESST 
collaboration demonstrated also that the pulse 
shape discrimination (PSD) technique can dis- 
tinguish WIMP-induced nuclear recoils from 
those induced by backgrounds by means of in- 
serting a scintillating foil, which causes some 
additional scintillation light for events induced 
by a-decay of 210 Po and thus shifts the pulse 
shapes of these events faster than pulses induced 
by WIMP interactions in the crystal [g]. 1 

However, as the most important issue in 
all underground experiments, possible residue 
background events which pass all discrimina- 
tion criteria and then mix with other real WIMP- 
induced events in our data sets should also be 
considered. Therefore, as a more realistic study, 
we take into account small fractions of residue 
background events mixed in experimental data 
sets and want to study how well the model- 
independent methods could extract the input 
WIMP properties by using these "impure" data 
sets and how "dirty" these data sets could be to 
be still useful. 



In this article, I focus on two properties of 
halo WIMPs: their mass m x and one-dimensional 
velocity distribution function /i(v). More de- 
tailed discussions can be found in Refs. 



2. Effects of residue background events 

In our numerical simulations based on the 
Monte Carlo method, while the shifted Maxwellian 
velocity distribution [Q, |I|] with the standard 
values of the Sun's orbital velocity and the Earth's 
velocity in the Galactic frame: vo ~ 220 km/s 
and v e = 1.05 vo, and the Woods-Saxon form 
for the elastic nuclear form factor for the spin- 
independent (SI) WIMP-nucleus interaction [ |Io[ 
^] have been use for generating WIMP-induced 
signals, a target— dependent exponential form 
for residue background events has been intro- 
duced [fill: 
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(2.1) 



Here Q is the recoil energy, A is the atomic 
mass number of the target nucleus. The power 
index of A, 0.6, is an empirical constant, which 
has been chosen so that the exponential back- 
ground spectrum is somehow similar to, but 
still different from the expected recoil spectrum 
of the target nucleus (see Figs. [[]); otherwise, 
there is in practice no difference between the 
WIMP scattering and background spectra. Note 
that, the atomic mass number A has been used 
here just as the simplest, unique characteristic 



parameter in the analytic form ( |2.1| ) for defin- 
ing the residue background spectrum for dif- 
ferent target nuclei. It does not mean that the 
(superposition of the real) background spectra 
would depend simply /primarily on A or on the 
mass of the target nucleus, m^. 

Note also that, firstly, the exponential form 



1 More details about background discrimination tech- 
niques and status see also e.g., Refs. [ph. 



(2.1) for residue background spectrum is rather 
naive; however, since we consider here only a 
few (tens) residue background events induced 
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by perhaps two or more different sources, pass 
all discrimination criteria, and then mix with 
other WIMP-induced events in our data sets of 
^(100) total events, exact forms of different 
background spectra are actually not very im- 
portant and this exponential spectrum should 
practically not be unrealistic. Secondly, our 
model-independent data analysis procedures re- 
quires only measured recoil energies from one 
or more experimental data sets with different 
target nuclei [[T]-^]. Hence, for applying these 
methods to future real direct detection data, the 
prior knowledge about (different) background 
source(s) is not required at all. 

Moreover, the maximal cut-off of the ve- 
locity distribution function has been set as 
Vmax = 700 km/s. The experimental threshold 
energy has been assumed to be negligible and 
the maximal cut-off energy is set as 100 keV. 
The background window (the possible energy 
range in which residue background events can- 
not be ignored, compared to some other ranges) 
has been assumed to be the same as the ex- 
perimental possible energy range. Note here 
that the actual numbers of generated signal and 
background events in each simulated experi- 
ment are Poisson-distributed around their ex- 
pectation values independently, and the total 
event number in one experiment is then the sum 
of these two numbers; both generated signal 
and background events are treated as WIMP 
signals in our analyses. Additionally, we as- 
sumed that all experimental systematic uncer- 
tainties as well as the uncertainty on the mea- 
surement of the recoil energy could be ignored. 

2.1 On the measured recoil spectrum 

In Figs. [I] I show measured energy spec- 
tra (solid red histograms) for a 76 Ge target with 
three different WIMP masses: 25 (top), 50 (mid- 
dle), and 250 (bottom) GeV. While the dotted 
blue curves show the elastic WIMP-nucleus scat- 
tering spectra, the dashed green curves indicate 
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Figure 1: Measured energy spectra (solid red 
histograms) for a 76 Ge target with three differ- 
ent WIMP masses: 25 (top), 50 (middle), and 
250 (bottom) GeV. The dotted blue curves are the 
elastic WIMP-nucleus scattering spectra, whereas 
the dashed green curves are the exponential back- 
ground spectra normalized to fit to the chosen back- 
ground ratio, which has been set as 20% here (plots 
fromRef. ftfo). 
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the exponential background spectrum given in 
Eq. ( |2.l[ ), which have been normalized so that 
the ratios of the areas under these background 
spectra to those under the (dotted blue) WIMP 
scattering spectra are equal to the background- 
signal ratio in the whole data sets. 5,000 exper- 1 
iments with 500 total events on average in each £ 
experiment have been simulated. 

It can be found here that, the shape of the 
WIMP scattering spectrum depends highly on 
the WIMP mass: for light WIMPs (m x < 50 
GeV), the recoil spectra drop sharply with in- 
creasing recoil energies, while for heavy WIMPs 
(m x > 100 GeV), the spectra become flatter. 
In contrast, the exponential background spectra 
shown here depend only on the target mass and 
are rather flatter/sharper for light/heavy WIMP 
masses compared to the WIMP scattering spec- 
tra. This means that, once input WIMPs are 
light/heavy, background events would contribute 
relatively more to high/low energy ranges, and, 
consequently, the measured energy spectra would 
mimic scattering spectra induced by heavier/ 
lighter WIMPs. 

2.2 On determining the WIMP mass 

Fig. U show the median values of the re- 
constructed WIMP mass and the lower and up- 
per bounds of the la statistical uncertainty by 
means of the model-independent procedure in- 
troduced in Refs. [Q] with mixed data sets from 
WIMP-induced and background events as func- 
tions of the input WIMP mass. As in Refs. [0], 
28 Si and 76 Ge have been chosen as two target 
nuclei. The background ratios shown here are 
no background (dashed green), 10% (long-dotted 
blue), 20% (solid red), and 40% (dash-dotted 
cyan) background events in the analyzed data 
sets. 2 x 5,000 experiments with 50 total events 
on average in each experiment have been sim- 
ulated. 

It can be seen here clearly that, since for 
light WIMP masses (m x < 100 GeV), due to 
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Figure 2: The reconstructed WIMP masses as func- 
tions of the input WIMP mass. 28 Si and 76 Ge have 
been chosen as two target nuclei. The background 
ratios shown here are no background (dashed 
green), 10% (long-dotted blue), 20% (solid red), 
and 40% (dash-dotted cyan) background events in 
the analyzed data sets. Each experiment contains 
50 total events on average. Other parameters are as 
in Figs. 1 (plot from Ref. [Q]). 

the relatively flatter background spectrum (com- 
pared to the scattering spectrum induced by WIMPs) 
or, in practice, some background sources in high 
energy ranges, the energy spectrum of all recorded 
events would mimic a scattering spectrum in- 
duced by WIMPs with a relatively heavier mass, 
the reconstructed WIMP masses as well as the 
statistical uncertainty intervals could be overes- 
timated. In contrast, for heavy WIMP masses 
(m x > 100 GeV), due to the relatively sharper 
background spectrum or e.g., some electronic 
noise, relatively more background events con- 
tribute to low energy ranges, the energy spec- 
trum of all recorded events would thus mimic 
a scattering spectrum induced by WIMPs with 
a relatively lighter mass. Hence, the recon- 
structed WIMP masses as well as the statis- 
tical uncertainty intervals could be underesti- 
mated. Nevertheless, Fig. || shows that, with 
~ 20% residue background events in the an- 
alyzed data sets of ~ 50 total events, the 1(7 
statistical uncertainty band can cover the true 
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WIMP mass pretty well; if WIMPs are light 
(m x < 200 GeV), the maximal acceptable frac- 
tion of residue background events could even 
be as large as ~ 40%. 

2.3 On reconstructing the one-dimensional : 
WIMP velocity distribution function 

In this section I show the median values 
of the reconstructed one-dimensional velocity 
distribution function of halo WIMPs 2 with its 
la statistical uncertainty by means of the model- 
independent method introduced in Ref . [|I|] with 
mixed data sets. As in Ref. [jlj], a 76 Ge nucleus 
has been chosen as our detector target for re- 
constructing /i(v); while a 28 Si target and a 
second 76 Ge target have been used for deter- 
mining m x . The background ratios shown here ; 
are no background (dashed green), 10% (long- 
dotted blue), and 20% (solid red) background 
events in the analyzed data set(s). (3 x) 5,000 
experiments with 500 total events on average in 
each experiment have been simulated. 

2.3.1 With a precisely known WIMP mass 

In Figs. |3| we first assume that the required 
WIMP mass for reconstructing f\ (v) has been 
known precisely with a negligible uncertainty. 
The horizontal bars show the sizes of the win- 
dows used for estimating /i(v) [jl|]. Five bins : 
have been used and up to three bins have been 
combined to a window, in order to collect more 
events in wider ranges and thereby reduce the 
statistical uncertainty on the reconstructed f\ (v). 

It can be seen that, since for heavy WIMP 
masses (m x > 100 GeV), the relatively sharper 
background spectrum contributes more events 
to low energy ranges, or, equivalently, to low 
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2 Note that, since the experimental maximal cut-off 
energy is fixed as 100 keV, for heavy input WIMP masses 
(m% > 250 GeV), one can reconstruct the velocity distri- 
bution function only in the velocity range v < 300 km/s. 

3 Note that, since the neighboring windows overlap, 
the estimates of f\ (v) at adjacent points are correlated. 



Figure 3: The reconstructed one-dimensional 
WIMP velocity distribution function for three dif- 
ferent WIMP masses: 25 (top), 50 (middle), and 
250 (bottom) GeV. The double-dotted black curves 
are the input shifted Maxwellian velocity distri- 
bution. The background ratios shown here are 
no background (dashed green), 10% (long-dotted 
blue), and 20% (solid red) background events in the 
analyzed data set. Each experiment contains 500 
total events on average. Other parameters are as in 
5 Figs. 1 (plots from Ref . [§]). 
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velocity ranges, the reconstructed velocity dis- , 
tribution would be shifted to lower velocities. 

0. 

For an input WIMP mass of 100 GeV and the 
background ratio of 10% (20%), the peak of 

the reconstructed velocity distribution function _ » 

l 

could be shifted by ~ 30 (60) km/s. In con- a , 
trast, for light WIMP masses (m x < 50 GeV), j 

o. 

the relatively flatter background spectrum con- 
tributes more events to high energy/velocity 
ranges, and the reconstructed velocity distri- ». 
bution would be shifted to higher velocities. 
Moreover, our simulation results indicate that, 
with an ~ 10% - 20% background ratio in the ■ 
analyzed data set of ~ 500 total events, one 
could in principle reconstruct the one-dimensional 

( 

velocity distribution function of halo WIMPs 

with an 6.5% (for a 25 GeV WIMP mass, ° 

20% background events) - ~ +38% (for a 250 1 . 
GeV WIMP mass, 10% background events) de- o 
viation. If the mass of halo WIMPs is ^(50 

( 

GeV), the maximal acceptable background ra- 
tio could even be as large as ~ 40% with a de- 
viation of only ~ +14%. 

2.3.2 With a reconstructed WIMP mass 

In Figs. |] the required WIMP mass for re- ». 
constructing f\ (v) has been reconstructed with 
other direct detection experiments. Note that, 

0. 

while the vertical bars show the la statistical j 
uncertainties on the reconstructed f\{v), taken s ' 
into account the statistical uncertainty on the ~ » 
reconstructed WIMP mass, the horizontal bars 
indicate here the 1 a statistical uncertainties on 

0. 

the estimates of the reconstructed v points, due 
to the uncertainty on the reconstructed m x . 

It has been found that, firstly, as shown in 
Figs. H for an input WIMP mass of 100 GeV, 
the reconstructed mass doesn't differ very much 
from the true value. Hence, the reconstructed 
fx (v) is approximately the same for both cases 
with the input and reconstructed WIMP masses. 
However, for light/heavy input masses (m x < / 
> 100 GeV), the reconstructed f x (v) with the 
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Figure 4: As in Figs. 3, except that the WIMP 
masses have been reconstructed by means of the 
procedure introduced in Refs. [^] (plots from 
Ref. §). 

reconstructed WIMP mass shift to relatively 
lower/higher velocities compared to the case 
with the input (true) WIMP mass. This effect 
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caused directly by the over-/underestimate of 
the reconstructed WIMP mass [^] : once the re- 
constructed WIMP mass is over-/underestimated 
from the real value, a transformation constant 
from Q to v will thus be under-/overestimated. 
Consequently, the reconstructed v points will 
be smaller/larger than the true values. 

In contrast to the first effect discussed in 
the previous subsection, this second effect draws 
the reconstructed WIMP velocity distribution 
(more strongly) to the opposite directions. Nev- 
ertheless, with an ~ 5% - 10% background 
ratio in the analyzed data sets of ~ 500 total 
events, one could in principle still reconstruct 
the one-dimensional WIMP velocity distribu- 
tion function with an ~ +7% (for 25 GeV WIMPs, 
10% backgrounds) - ~ +16% (for 250 GeV 
WIMPs, 5% backgrounds) deviation. If m x < 
100 GeV, the maximal acceptable background 
ratio could even be as large as ~ 20% with a 
deviation of only ~ +9%. 

3. Summary 

In this article we reexamine the data analy- 
sis methods introduced in Refs. [jlj ||] for deter- 
mining the mass and one-dimensional velocity 
distribution function of halo Dark Matter par- 
ticle from measured recoil energies of direct 
detection experiments directly, by taking into 
account small fractions of residue background 
events, which pass all discrimination criteria 
and then mix with other real WIMP-induced 
events in the analyzed data sets. 

Our simulations show that, with a back- 
ground ratio of ~ 5% - 10% in data sets of ~ 
500 total events, the one-dimensional WIMP 
velocity distribution function can in principle 
be reconstructed with an ~ —5% — h20% sys- 
tematic deviation; whereas with a background 
ratio of even ~ 10% - 20% in data sets of only 
~ 50 total events, the 1 a statistical uncertainty 



band of the reconstructed WIMP mass can cover 
the true value pretty well. 
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